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Abstract 


Cave geotechnical studies have been the key to meeting the requirements of Brazilian 
environmental legislation for the conservation of speleological heritage in mining 
areas. This paper presents a methodology that classifies iron caves according to their 
susceptibility to structural instability called the Cave Geomechanical Index (CGI). 
This index combines four variables: (1) Rock Mass Rating (RMR), Bieniawski’ s 
geomechanical variable, which classifies the quality of the rock mass hosting the 
cave; (2) Hydraulic Radius (HR), an engineering variable that allows the dimension 
of the span to be evaluated; (3) Ceiling Shape (CS), a speleological variable that 
indicates whether the ceiling geometry of the cave spans is favorable or unfavorable 
for block formation, and (4) Ceiling Thickness (CT), a geotechnical variable that 
represents the depth between the ceiling of the cave and the surface of the ground 


Corresponding Author: regarding auto-support issues. The CGI was developed, applied and calibrated 

Marcelo Roberto Barbosa over four years, by monitoring 63 spans from 27 caves adjacent to active iron ore 

Vale S.A., Nova Lima, Minas Gerais, Brazil. mines in Carajas, ParA state, Brazil, that had been authorized to be eliminated. 
Email:marcelo.roberto.barbosa@vale.com This geomechanical classification system proved to be easy to implement and its 


results showed that 76% of the spans with breakdown occurrences in the caves were 
classified as high or very high susceptibility to structural instability, while 94% of 
the spans classified as low susceptibility did not show any signs of physical damage. 
The CGI is discussed with the focus on improving stability studies, predictability 
of cave breakdown mechanisms and geotechnical risk analysis of iron caves near 
mining operations. 


Keywords: Iron cave, Stability assessment, Cave breakdown, Cave geomechanics, Risk 
analysis, Carajas 


Introduction cal-scientific studies are required in the environmen- 


In Brazil, natural underground caves are protected by tal licensing processes that regulate the preservation 


federal environmental legislation and several techni- _ of Brazilian speleological heritage. Caves abound in 


iron formation terrains so in recent years researchers 
have intensified a balance between the iron mining 
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Most of the iron caves in Brazil are found in two 








large ferruginous geosystems: the Iron Quadrangle 
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region of Minas Gerais State, southeastern Bra- 
zil, and the Carajas ridge, of Para State, in north- 
em Brazil (Valentim & Olivito 2011). Currently, 
the National Speleological Database of the federal 
environmental agency, Chico Mendes Institute of 
Biodiversity ICMBio), has about 4,300 iron caves 
registered, most of which are located in the Carajas 
ridge of Para State. This is due to a greater effort of 
cave prospecting in areas where iron ore extraction 
is higher, demonstrating the strong link between the 
evolution of speleological knowledge and the iron 
ore mining industry. 


Federal Decree 6.640/08 required that caves must be 
classified according to 11 physical, biological, ecolog- 
ical, and/or historical-cultural attributes, to establish 
their level of significance (importance for the spele- 
ological ecosystem). About these attributes, scientific 
studies are carried out to advance the knowledge on 
each of the themes allowing the correlation and the 
precise classification of the significance. In the archae- 
ological context, extremely important in the recovery 
of the ancient history of past populations, all sites are 
duly rescued for Brazilian museums. 


The mentioned Decree also regulates the different 
ways to license cave impacts, depending on their lev- 
el of significance, or requiring studies that guarantee 
that any anthropic actions will not compromise the 
physical integrity of the caves (Brasil 2008). 


During these often long-term studies within the licens- 
ing issues, another regulation, CONAMA 347, still 
requires that a buffer zone of 250 m should be pro- 
tected around each cave (MMA 2004). In the context 
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of mineral industry, this regulation directly interferes 
with mining operations, blocking large areas with the 
substantial amount of mineral reserves, which may 
lead to the stopping of mines and/or making the mine 
design more complex increasing production costs. 


Therefore, it is important and urgent for the mining 
companies to understand the geomechanical mech- 
anisms for the structural stability of caves in mining 
regions, which will allow greater security for the op- 
erations and leading confidence for the environmental 
agencies that may speed up the licensing of the areas. 


This paper introduces the Cave Geomechanical 
Index (CGI), a novel system for classifying the 
susceptibility of iron cave spans to structural insta- 
bility. This index was based, conceptually, on the 
Bieniawski (1989) rock mass geomechanical quali- 
ty rating (RMR) system, but is focused exclusively 
on natural iron caves, considering other essential 
variables for stability analysis. 


A simplified development flow of the main activi- 
ties for the development of CGI is shown in Fig. 1. 


The CGI was applied to 63 spans of 27 strategically 
selected caves, being the most representative in spe- 
leological terms and the closest to the mining areas. 
They occurred near the N4 and N5 mines in Carajas 
ridge which were used for development, adjustments 
and calibrations of the index between 2015 and 2018. 


Iron Caves 
Iron caves are poorly studied when compared with 
caves in carbonate rocks or other types of terrain. 
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Figure |. Simplified development flow of the main activities of the Cave Geomechanical Index (CGI). 
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Generally, iron caves are shallow with thin ceilings 
(20 m average), with their development coincident 
with the slope of the overlying land surface. The 
spans of iron caves are small (20 m length x 10 m 
width average), and irregular in shape, with struc- 
tural control and discontinuities observed on their 
walls and ceilings (Pild & Auler 2009; Pild et al. 
2015). White & White (1969) and White (2012) re- 
ported that due to the natural evolution of caves and 
the breakdown along discontinuities, it is common 
to observe fragments and blocks of rocky material 
on the cave floors. 


The genesis of iron caves is attributed to four main 
processes of generation of empty spaces in rock 
masses: erosion, leaching, dissolution, and biogen- 
esis (Simmons 1963; Vann 1963; Moss 1965; Mor- 
ris 1985; Pinheiro & Maurity 1988; Buchmann et 
al. 2009; Calux 2013; Dutra 2013). Studies of these 
processes are keys to understanding the geotechni- 
cal issues and stability analysis of this type of cave. 


As an example, Fig. 2 shows the topographic map 
and photos of cave N4E 0026, one of the studied 
caves. 
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Geotechnics Applied to Speleology 

The vast majority of cave stability studies reported 
in the literature have been undertaken in carbonate 
terrains (Waltham 2002; Parise & Trisciuzzi 2007; 
Hatzor et al. 2010; Szunyogh 2010; Gutiérrez et 
al. 2014; Parise et al. 2015; Jorda-Bordehore et al. 
2016, 2017; Andriani & Parise 2017; Fiore et al. 
2018). Carter and Miller (1995) observed that in 
carbonate rocks there is some relationship between 
ceiling width and cave stability. Criteria based on 
RMR (Bieniawski 1989) have also been used in the 
geotechnical mapping of these areas (Siegel & Mc- 
Crackin 2001). 


Some ideas about structural instability from car- 
bonate caves however can be transferred to the iron 
caves environment, leading to recent studies aim- 
ing to predict irreversible negative impacts on spe- 
leological heritage on mining sites (Sanchez 2007; 
Brandi 2015a, b; Araujo et al. 2016; Rend 2016; 
Valentim 2016; Dutra 2017; Lacerda 2017; Santos 
Junior 2017; Novas et al. 2017, Brandi 2019a). 


Also regarding the ferriferous cave environment, 
Reno (2016) developed a cave zoning system to 
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Figure 2. Photos of cave N4E_0026 (one of the studied caves). 
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support geotechnical studies. Using this approach, 
the geotechnical analysis of the spans of some iron 
caves was undertaken based on traditional geome- 
chanical classification systems adapted to caves. 
In line with this research, the works of Araujo et 
al.(2016) and Valentim & Olivito (2016) corroborat- 
ed the evolution of the geotechnical theme applied 
to speleology in iron caves from, respectively, the 
survey of geomechanical parameters with a 3D laser 
scanner, and the proposed classification of a cave 
through the method of geomechanical classifica- 
tion of rock masses adapted to caves. Dutra (2017) 
evaluated the geotechnical susceptibility of two iron 
caves in the Serra do Gandarela region by mapping 
parameters as geomechanical classification of rock 
mass, structural discontinuities, ceiling thickness, 
water input points (drips, speleothems, basins), sed- 
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iment cones, root incidence and characterizing the 
rainfall and water infiltration into the caves. 


Study Area 

The study area is located at the N4 and N5 iron 
mines in Carajas ridge southeastern of Para State. 
The selection of the caves took into account their 
physical representativeness of the Mineral Province 
of Carajas as horizontal projection (from 5 to 200 
m) and volume (from 5 to 1000 m°), and their prox- 
imity to active mining areas (Mines N4EN, N4WS 
and NSS). Twenty seven caves were selected which 
have been monitored over 4 years between 2015 
and 2018 (Fig. 3). 


Table 1 shows the level of significance, lithology 
and morphometric data of the 27 caves of this study. 
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Figure 3. Study area at the iron mines N4EN, N4WS and NSS in the Carajas ridge, Para State. Note the 
location of the 27 caves (yellow bullets) selected for this study and their proximity to the mine opera- 
tions (Worldview image (2019); UTM coordinates). 





Volume 3 / Issue 2/ pages(134-161) e-ISSN: 2588-7343 





p-ISSN: 2645-4661 137 








Table 1. Level of significance, lithology and speleometry data of the 27 caves of the study (BIF = Banded 


Iron Formation, Canga = Ferruginous Crust). 




















































































































Cave Eve ot vent Lithology Plan Lenght (m) Sees Area (m?) Volume (m7?) 
cance (m) 
N4E 0015 High BIF and Canga 65 6 M33) 547 
N4E_0016 High BIF 30 32) 113 114 
N4E_ 0021 High BIF 41 6.6 164 330 
N4E 0022 Maximum BIF and Canga 175 12.4 1082 750 
N4E 0023 Maximum BIF and Canga 65 5.6 274 POP 
N4E 0026 High BIF and Canga 162 18 451 695 
N4E_0027 High BIF and Canga 14 35) 24 29 
N4E_0028 High BIF 25 D5) 81 113 
N4E_0037 High BIF and Canga ne) 2) 19 18 
N4E_0039 High BIF 30.6 3.6 69 67 
N4E_0044 High BIF and Canga Poul 1.1 48 331 
N4E_0061 High BIF and Canga 335 1 112 134 
N4E 0086 High BIF and Canga USS) ei 32) 3h) 
N4WS _ 0001 High BIF 13 44 64 68 
N4WS_0008 High BIF and Canga 79 Q 445 876 
N4WS_0022 High BIF 10.5 2) Sil 91 
N4WS_0050/51 High BIF and Canga MDS 4.8 178 207 
N4WS_0059 High BIF and Canga 16 Be 52 67 
N4WS_0060 High BIF 12.5) 7 15 8 
N4WS_0065 High BIF 13355) 1.8 15 12 
N5S_0001 High BIF and Canga 16 3 31 40 
N5S_0011 High BIF and Canga 104 11 471 898 
N5S_0023 High BIF 26.5 2S) 32 19 
N5S_0071 High BIF and Canga 8.3 1.6 17 19 
N5S_0072 High BIF and Canga a5 2) val 82 
N5S_0073 High BE aa 20 0.6 43 28 
= Canga 
N5S_0075 High BIF and Canga 32 3 76 59 








Geology 

The geology of the area is characterized by iron 
deposits distributed along a structure known as the 
syncline or Carajas fold, which is grouped into sets 
called Serra Norte, Serra Sul, Serra Leste, and Sao 
Félix. More specifically, the caves considered in this 
study are developed in iron lithologies of the Cara- 
jas Formation, which is part of a Carajas Domain’s 
greenstone belt. 


The stratigraphic unit includes Precambrian banded 
iron formation (BIF) of oxide facies with jaspilites 
(Lateritic Cover and Carajas Fm) interdigitated by 
moderately to severely deformed metabasalt and 
acid metavolcanic packages (Parauapebas Fm), me- 
tavolcanic sedimentary rocks into green to amphib- 
olite sub-schist facies (Igarapé Pojuca Fm), and do- 
main of metasedimentary rocks (Aguas Claras Fm), 
complementing with the Proterozoic granites (Po- 
juca, Igarapé Gelado and Central), (Justo & Lopes 
2014) (Fig. 4). 


The high iron content of Carajas Formation lith- 


ologies enables the formation of ferruginous crust 
commonly called “canga”. This Cenozoic deposit 
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occurs in a residual form when formed in situ by the 
leaching of jaspilite chert bands or other minerals/ 
constituents of other rich ferruginous rocks; or set- 
tles in detrital form when the source rock (rich in 
iron) is disaggregated, transported, and subsequent- 
ly cemented. These deposits are up to 20 m thick 
and regularly hard, and provide excellent support 
for underground cave spans. 


Cave Geomechanical Index (CGI) 

The development of indexes that evaluate the geo- 
hazard impacts are not new for carbonate karst 
environments. Van Beynen & Townsend (2005) 
initially proposed the “Karst Disturbance Index”, 
later revised by North et al. (2009). Beynen et al. 
(2012) created a more widely index called “Karst 
Sustainability Index” the standardizing metric of 
sustainable development practices in karst settings 
that takes into account 25 indicators related to the 
environment, economy, and society. A high-level 
review was also provided in Gutiérrez et al. (2014) 
reviewing the main natural and human-induced haz- 
ards and their indexes characteristic of carbonate 
karst environments, focusing on sinkholes, floods 
and slope movements. It was quite relevant this lat- 
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Figure 4. Regional geology of the mineral province of Carajas with the approximate location of 
study caves (adapted from Dardenne and Schobbenhaus, (2001)). 
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ter work highlighting 17 main causes of degrada- 
tion and a classification for their effects, raising 3 of 
them related to engineering and mining operations. 


The differences and peculiarities between caves in 
carbonate karst and ferruginous karst environments 
do not allow the direct use of all the results and tools 
produced for the carbonate environments, although 
many of the studies in limestone caves, especially 
those dealing with structural stability, have been 
used and cited for the present study. 


The Cave Geomechanical Index (CGI) is a geome- 
chanical classification system specifically for iron 
ore caves near mining areas that are required to 
have their physical integrity preserved accordingly 
to current Brazilian legislation. 


Definition of Formula and Variables 

The CGI development model was inspired by Bi- 
eniawski’s (1989) geomechanical classification 
methodology, with the sum of weighted variable 
classes according to ranges or interval of analysis 
and whose result indicated the geomechanical qual- 
ity of the evaluated material. The variables selected 
for the CGI and its classes cover several disciplines 
and were based on the experience and knowledge of 
the authors in geostructural mapping and monitoring 
since 2014 in caves located near active mining areas 
in Carajas, Para State and in the Iron Quadrangle re- 
gion, Minas Gerais State, both in Brazil, as well as 
based on the scientific literature. Three recent articles 
are representative of this gain of experience in mon- 
itoring iron caves in near mine operations: Barbosa 


(2018), Barbosa et al. (2019) and Brandi (2019b). 


Four variables were considered: (1) Rock Mass Rat- 
ing (RMR), from Bieniawski (1989), a geomechan- 
ical variable that classifies the quality of the rock 
mass in which each cave is hosted; (2) Hydraulic 
Radius (HR), an engineering variable that considers 
the ratio between the area and the perimeter of a 
typical section of a span; (3) Ceiling Shape (CS), a 
variable that indicates whether the cave ceiling ge- 
ometry is favorable or unfavorable to form blocks; 
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and (4) Ceiling Thickness (CT), a geotechnical vari- 
able that characterizes the depth of the cave span 
ceiling concerning the ground surface. 


The system is governed by an equation that consid- 
ers the sum of weighted variable classes which leads 
to a classification of the geomechanical quality of 
the caves’ spams, measuring their susceptibility to 
structural instability. 


The variables RMR, HR and CT (quantitative), 
and CS (qualitative), induce different effects on 
the classification of the CGI, depending on the 
combination of these variables in its equation. Each 
one contributes to different impacts on the results 
and so on its classification. The higher the RMR 
and CT values, each of them contributes to a better 
classification of the GGI, taking the classification 
to very low susceptibility to structural instability. 
About HR, the opposite is true, the higher its value, 
it will contribute to a worse classification of the 
CGI, leading to very high susceptibility to structural 
instability. Regarding CS, the CGI classification 
gets worse or improve according to the type of 
ceiling shape, more or less unstable. 


The CGI formula is demonstrated in Equation 1: 
CGI= a RMR + BHR +y CS +6 CT ......000 Equation 1| 
Where: 

a RMR = Weight assigned to Rock Mass Rating 
(Bieniawski, 1989) 

B HR = Weight assigned to the Hydraulic Radius 
(Laubscher, 1990) 

y CS = Weight assigned to Ceiling Shape 

5 CT = Weight assigned to Ceiling Thickness 


The sum of the weights of the 4 variables was de- 
fined as the maximum value 100 (one hundred), 
which represents the best scenario of the very low 
susceptibility to structural instability, whereas, on 
the other hand, the worst scenario is the minimum 
value 0 (zero), indicating very high susceptibility 
to structural instability. Table 2 shows the five CGI 
classes and their respective levels of susceptibility 
to structural instability of the spans. 
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Table 2. Five CGI classes and their respective levels 
of susceptibility to structural instability of the spans 
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Variable classes 
Based on the acquired knowledge, mainly in the 
field observation of the caves immediately af- 





























Susceptibility to Structural Insta- fe ter the detonations, it was possible to recognize 
bility of Cave Span and propose the best and worst classes for each 
Very Low CGI > 80 variable according to Table 3, which will be de- 
Eon EUS CoH Sy scribed in this chapter, always exemplifying its 
Moderate SS EG relations with cave N4E 0026, one of the 27 
ee BU SCC Sst studied caves. 
Table 3. CGI variable classes 
Rock Mass Rating (RMR) Hydraulic Radius (HR) Ceiling Shape (CS) Ceiling Thickness (CT) 
¢ Very Good ¢ Small ¢ Arch ¢ Large 
* Good ¢ Regular ¢ Planar ¢ Regular 
¢ Regular ° Large ¢ Inverted Arch ¢ Small 
* Poor 
¢ Very Poor 
Legend: Best Class Regular Class Worst Class 











Rock Mass Rating (RMR) 

Since the 1940s, several authors have designed geo- 
mechanical indexes to predict the behavior of rock 
masses, particularly when exposed to anthropic ac- 
tions such as mining and tunneling (Edelbro 2006). 


For the composition of the CGI, 3 traditional sys- 
tems were considered for their approaches to mass 
quality and application in underground environ- 
ments, as is the case of the shallow iron caves of the 
study area. Barton’s Q (1974), Bieniawski’s RMR 


N (RMR) BIENIAWISKI 


(Q) BARTON 


ol Poor 





BB Poor Rock 





(1989) and Laubscher’s MRMR Systems (1990) 
were tested. The tests were performed in the 27 
study caves, making a correlation with the geostruc- 
tural field mapping. The results were very similar, 
but the Bieniawski RMR System was perceived as 
the most adherent (De Paula et al. 2018), and the 
one that provided better geomechanical contact de- 
tails to compose the CGI. 


Fig. 5 shows an example of the mapping result in 
the 3 systems for cave N4E _ 0026. 


Figure 5. Geomechani- 


cal classification maps 
of cave N4E_ 0026, fol- 
lowing Bieniawski RMR 
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(1989), Barton Q Sys- 
tems (1974) and Laub- 
scher (1990). 
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Bieniawski proposed the RMR system in 1973, and 
since then it has been refined until its last revision 
in 1989. It is based on the analysis of six parameters 
(uniaxial compressive strength of intact rock, rock 
quality designation, joint or discontinuity spacing, 
joint condition, groundwater condition and joint ori- 
entation) related to the rock mass. Each parameter 
is given a weight and the sum of the weights is the 
result of the classification system. Before being ap- 
plied, the rock mass must be divided into structural 
sectors with similar characteristics, and the classifi- 
cation applies to each of these sectors. 
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After obtaining the weight values for the parame- 
ters related to the geomechanical quality of the rock 
mass and the discontinuities, Bieniawski suggested 
performing a correction related to the orientation of 
the discontinuities according to their influence on 
the structure stability. The result of this geomechan- 
ical quality classification system of the rock mass 
was divided into five classes with specific intervals 
and was used for the composition of the CGI, as 
follows: Very Good (100 to 81), Good (80 to 61), 
Regular (60 to 41), Poor (40 to 21), and Very Poor 
(20 to 0), (Table 4). 


Table 4. Rock mass rating (RMR) from Bieniawski (1989) used in the CGI. 























Sum 100 - 81 80 - 61 60 - 41 40-21 <20 
Class number I IL Ill IV Vv 
Description Very Good Good Regular Poor Very Poor 














Hydraulic Radius (HR) 

Several researchers have proposed indexes/param- 
eters about horizontal sections when assessing the 
stability of underground spaces. Two main con- 
cepts were: ‘’*Span Width’, developed by Barton 
et al. (1974) and Hutchinson & Diederichs (1996), 
and ‘’Hydraulic Radius’, developed by Laubscher 
(1990). 


Span Width 

Milne (1997) reports that the span width support 
condition has been employed for decades and one 
of the earliest records goes back to Terzaghi (1949). 
This method is mainly used in tunnel designs and 
can be effectively related to stability since the width 
of the opening defines the distance between the 
main supporting elements of a tunnel, which are the 
sidewalls. 


Hydraulic Radius 

Hydraulic radius has been used in the assessment of 
the stability of underground structures since 1977 
and was first applied by D. H. Laubscher (Milne 
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1997). Milne points out that Hydraulic radius was 
initially used in fluid dynamics to compare the flow 
in square section tubes with flow in circular section 
tubes. The results explained why a parameter ini- 
tially applied in hydraulics is related to the stability 
of underground excavations. 


The Hydraulic radius is obtained by the ratio be- 
tween the area and the perimeter of the analyzed 
section. This ratio can also express the distance 
between the center of the underground space and 
the supporting structures, as well as when dealing 
with irregularly shaped excavations, and when the 
structure length is up to three times the width since 
this parameter considers all faces of the excavation 
as capable of providing effective support (Milne 
1997). 


According to Laubscher (1990), the relationship 
between the hydraulic radius and the geome- 
chanical classification index (Mining Rock Mass 
Rating — MRMR) was proposed for the analysis 
of self-portable of spans in underground mines. 
Based on data from 11 mines in four countries, an 
abacus was generated, classifying points as stable, 
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transitional or unstable. Jorda-Bordehore (2017) 
evaluated the stability of 137 wide-spanned nat- 
ural caves based on the empirical approach, based 
on Barton’s Q system (1974), and span width, pro- 
posing a stability curve. 


The delimitation of the “effective span” is based on 
large areas (> 4 m’) with high ceilings (> 1,5 m) 
trying to reach the entire cave, however, some parts 
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of the cave may not be considered due to small areas 
end the low height of their ceilings. 

For the CGI composition, this geometric variable will 
be included according to 3 classes (Small, Regular and 
Large), which will represent their respective sizes es- 
tablished through statistical study in the metric unit. 


A map was prepared for each of the 27 studied caves 
with the division of spans considered in the calcula- 
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Figure 6. Hydraulic radius (HR) map of cave N4E_ 0026. Note that each cave span is identified 


with a sequential number. 


tion of the respective hydraulic radius. As an exam- 
ple, a span map of cave N4E_0026 is shown in Fig. 6. 


Ceiling Shape (CS) 

Peck et al. (2013) stated that the shape of the ceiling 
is a key factor in assessing the stability of under- 
ground excavation, along with the geological char- 
acteristics and general condition of the ground. Peck 


Volume 3/ Issue 2/ pages(134-161) e-ISSN: 2588-7343 p-ISSN: 2645-4661 y 


et al. (2013) compared the performance of straight 
(flat) and arched (arch) shapes in shallow excava- 
tions, concluding that areas of instability are usually 
associated with flat geometries and that stable ceil- 
ing shapes are similar to the post-rupture configura- 
tion -i.e. arches. Parise and Lollino (2011) evaluated 
the geological failure mechanisms in an abandoned 
calcarenite mine and recorded the results of the 


5) 





collapse, noting that the failure points have arched 
shapes. Diederichs & Kaiser (1999) suggested that 
discontinuities and intersections reduce the rock’s 
ability to form a stable arch, extending the unsta- 
ble zone. Similarly, Fiore et al. (2018) evaluated the 
collapse susceptibility of underground calcarenite 
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mines in Italy and recognized that ceiling geometry 
was an important factor, being the inverted arch the 
worst scenario. 


To determine the ceiling shape for the CGI, each of 
the 63 spans of the 27 study caves was analyzed in 


INVERTED ARCH 











their ceiling shapes in 3D laser sections, confirmed 
by field surveys and speleometric map profiles. 
Three classes were considered (Arch, Planar and 
Inverted Arch) (Fig. 7). 


Ceiling Thickness (CT) 

Several researchers have suggested that, in the case 
of shallow excavations, the closer to the surface, the 
less stability these structures present (Carter 2014). 
Yardimci et al. (2016) reported an attempt to evaluate 
stability using the Crown Pillar method. This meth- 
od establishes a relationship between geomechanical 
classification systems, such as RMR (Bieniawski 
1989) and Q (Barton et al. 1974), and the geome- 
try of underground spaces, considering their ceiling 
thickness. Carter & Miller (1995) proposed accept- 
able risk exposure guidelines, dividing underground 
mine sections into quality classes. This is a practical 
method that relies on back-analysis of historical fail- 
ure cases. Carter (2014) developed an equation that 
correlates geometric factors and ceiling thickness to 
calculate the scaled span (Cs) factor. Thus, small ceil- 
ing thicknesses generate high Cs values, which im- 
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Figure 7. Type sections of the 3 classes 
of the Ceiling Shape variable, extracted 
from the 3D laser sections of the study 


Caves. 


plies instability. Parise & Lollino (2011) relate ceil- 
ing thickness to the stability of underground spaces. 
Fiore et al. (2018) studied shallow structures stating 
that ceiling thickness is fundamental for the evalua- 
tion of underground space stability. 


To determine the ceiling thickness for the CGI, la- 
ser scanning was used to acquire topographic data 
with LIDAR technology producing 3D maps (inside 
and outside the caves). From the point cloud, a cave 
ceiling elevation model was generated and subtract- 
ed from the digital terrain elevation model to obtain 
the ceiling thickness values. 


Through the point-to-point statistical analysis, it 
was possible to calculate the average value for ceil- 
ing thickness of each span. Fig. 8 shows the ceiling 
thickness map for cave N4E_0026 in a 3D model. 


Value Ranges of Variable Classes (statistical 
study) 

The class value ranges of the 4 CGI variables were 
obtained as follows. Among the 3 quantitative vari- 
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Ceiling Thickness (m) 
0 4 8 12 14| ~~~ Span axis 
SS 9} E Cave entrance 


Figure 8. Ceiling thickness map for cave N4E_0026 in a 3D model. 


ables, RMR used the Bieniawski (1989) values, and COMMON ceiling shapes in the region. The qualifica- 
HR and CT used a statistical study developed in tion of the classes (best, regular, worst) is based on 
the MiniTab software from geostructural mappings the experience of the technical project team since 
data collected in the 63 spans of the 27 study caves 2015 (expert opinions). 

(Pareto distribution diagram with cumulative rela- 

Table 5 shows the range results for each of the 4 


tive frequency). The qualitative variable (CS) was 
CGI classes. 


obtained from the experience of mapping the most 
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Table 5. Class value ranges for the 4 CGI variables. 








Rock Mass Rating Hydraulic Radius Ceiling Shape Ceiling Thickness 
(RMR) (HR) (CS) (CT) 
¢ Very Good (81 - 100) ¢ Small (0.0-0.91m) ¢ Arch ¢ Large (> 7.65 m) 

















Legend: Best Class 

















Table 6. Data collected from the 4 CGI variables for the 63 spans of the 27 study caves. 































































































A ae 4 A Ceiling Thick- 
Gave Horizontal Projection Span | RMR Hydraulic Radius Ceiling Shape fies 
(m) (m) 
(m) 
1 50 1.78 Inverted Arch 6.62 
N4E_0015 65 
2} 12 2.85 Arch 4.59 
i 61 1.95 Planar 7.76 
N4E_0016 30 
2 62 1.08 Planar 10.51 
1 66 152) Planar 6.22 
N4E_0021 Al 
2 39 2.44 Planar 8.32 
1 ee E38) Planar 6.4 
2, 50 SP) Planar 6.16 
3 61 1.04 Arch 4.68 
N4E_0022 175 
4 37 2335) Inverted Arch 7.97 
5 34 1.11 Arch 8.40 
6 48 2.42 Inverted Arch 5.64 
1 44 2.38 Inverted Arch O71 
N4E_0023 65 
2; 39 1.82 Arch 10.39 
if 29 1.46 Arch lil 
2 42 Lay Planar 6.54 
3 64 0.79 Arch 5.20 
N4E_0026 162 
4 36 1.85 Planar 6.31 
5) 32) DSI Arch 8.89 
6 25) PB Inverted Arch 10.67 
N4E 0027 14 1 66 0.90 Planar 4.82 
N4E 0028 25 il 44 1.54 Inverted Arch 3.00 
N4E 0037 il? 1 78 0.79 Arch 2.70 
1 58 0,62 Arch 4.73 
N4E_ 0039 30,6 2 68 0.96 Inverted Arch 3.51 
3 68 1.40 Arch 3.79 
































































































































Gave Horizontal Projection Span | RMR Hydraulic Radius Ceiling Shape avon 

(m) (m) ae 
1 68 1.26 Inverted Arch 2.05 

N4E 0044 Zl 
a Sw 0.76 Planar 2.96 
N4E 0061 35 1 2il 1.87 Arch 3.50 
N4E 0086 155) 1 22 Nee Arch 3.74 
N4WS 0001 13 1 66 1.65 Inverted Arch 6.77 
1 65 1.35 Planar 7.64 
2) 30 2.60 Arch 8.56 
N4WS_0008 79 3 41 0.97 Planar 11.07 
4 30 2.44 Arch 7.64 
5 28 1.64 Planar 8.77 
N4WS_0022 10,5 1 73 0.66 Planar 2.00 
1 1S 0.92 Planar 2.96 
2 68 1.96 Inverted Arch 4.26 
N4WS_ 0050/51 Hydra 20,5 3 64 1.09 Planar 5.49 
4 83 0.73 Planar 6.66 
5 64 0.38 Inverted Arch 6.93 
N4WS_0059 16 1 65 0.95 Planar leafy 
N4WS_0060 IPE) 1 53 0.81 Planar 1.54 
N4WS_0065 13.5) 1 Sul 0.78 Arch 5.48 
1 61 0.87 Planar 3.40 
N5S_0001 16 2 71 0.43 Arch 4.33 
3 55) 0.6 Arch 3.31 
75 DS Planar 2.03 
2 1S 327] Arch 4.21 
3 60 2? Arch 7.28 
4 42 eS Planar 7.02 

N5S_0011 104 
5 59 0.91 Planar 8.17 
6 58 1.74 Planar 8.89 
7 42 0.86 Planar 9.36 
8 63 0.91 Planar 2.47 
1 69 0.42 Arch 3.18 

N5S_0023 26,5 
2 64 0.5, Planar 3.50 
NS5S_0071 8,3 1 53} 0.92 Arch 0.96 
NS5S_0072 2s) 1 3) 1.69 Inverted Arch 2.54 
NS5S_0073 20 1 3) 1.28 Planar 2M 
1 40 1.45 Planar Sil 

NS5S_0075 32 
2 55) sil Planar 4.24 





























Table 6 shows the result of the mapping of the  pyg Se ht: nC oe | 


CGI variables in the 63 spans of the 27 study The weights oe ai osauble were etre based on 


caves. the experience of the team of geologists and geotech- 


nicians who followed the geostructural mapping and 
monitoring project since 2015 (expert opinions), in- 
cluding the caves that could be monitored until their 
collapses. This situation is the most favorable as it 
allows technicians to make the closest correlations to 
reality between the approach of the controlled mine 
until the total cave collapse. This makes the proposed 
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weight scenarios much more assertive. To choose the 
best distribution of the proposed weights for each 
CGI variable, four possible scenarios were thought 
out in order to allow simulations with the CGI for- 
mula and thereby check the resulting susceptibility 
classifications that came closest to the real situation 
based on field experience. 


Table 7. Four scenarios with the proposed weights for simulation. 




































































Scenario 1 Scenario 2 Scenario 3 Scenario 4 
CGI Variable Class (proposed wei- (proposed wei- (proposed wei- (proposed wei- 

ght) ght) ght) ght) 
Very poor 0 0 0 0 

; Poor 6.25 12.5 21 21 

Rock Mass Rating 
Regular 12.5 25 39 39 
(RMR) zi 7 2 z 

700d 18.75 Ses 54 54 

Very good 25 50 60 60 

Small 25 28 20 25 

Hydraulic Radius (RH) Regular 12.5 14 12 15 
Large 0 0 0 
a ‘ Small 0 0 0 

Ceiling Thickness 

on Regular 12.5 ZS 3 %) 
oS) Large 25 5 5 5 

Arch 25 17 15 10 
Ceiling Shape (CS) Planar 1255 8.5 9 6 
Inverted Arch 0 0 0 0 

Legend: Best Class Regular Class Worst Class 








Table 7 shows the four scenarios with the distribu- 
tion of proposed weights for all classes of CGI vari- 
ables for simulation. 


An analysis of Table 7 showed that, excluding Sce- 
nario | of same weights for all best classes of vari- 
ables, it is noted that the variable Rock Mass Rat- 
ing (RMR) has a higher weight value over the other 
variables since it is the only one that reveals the 
geomechanical quality of the rock mass of the span, 
including the structural and hydrological features. 
The other CGI variables are equally important, but 
with lower weight values. 


Theoretical Weights , Better Adherence to Laub- 
scher’s Stability/Instability Diagram 
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The CGI system was simulated with the proposed 
weights of the four scenarios for the 63 spans of 
the 27 study caves and the results were inserted in 
Laubscher’s stability/instability diagram (Laubscher 
1990). This diagram is a world-renowned geome- 
chanical tool divided into three curves of stability 
classifications: unstable, transition and stable, which 
empirically assess the stability of underground areas. 
This fact sought to associate the CGI with mature 
tools widely applied and cited in the literature. 


All CGI classification results inserted into the Laub- 
scher’s diagram allowed to infer other curves result- 


ing in five curves of instability, creating a combined 
diagram called “Laubscher adapted with CGI’, with 
CGI ratings of spans of very high susceptibility to 
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Laubsher 


Laubsher Limit 
© Stable 
A Transition 


gw Unstable 


cGl 


Study Cave (with CGl) 


= NAWS_0001 
= NAWS_0008 

~ NAWS_0022 

= NAWS_0050/51 
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* NAWS_0065 
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Figure 9. Laubscher’s stability/instability diagram adapted with CGI. Note the five curves of susceptibility to instability 


inferred to CGI ratings. 


structural instability were within the unstable Laub- 
scher’s zone; the spans with very high to high sus- 
ceptibility CGI ratings were within the Laubscher’s 
transition zone, and finally the moderate, low, and 
very low CGI ratings of the spans were within the 
stable Laubscher’s zone (Fig. 9). 

This check of the proposed weights adherence 


through the Laubscher’s adapted diagram is part of 
the CGI development flow for the definition of the 
theoretical weights to be effectively used. Of the 
four simulated scenarios, scenario 4 was the most 
adherent to the adapted diagram, resulting in that 
of the 63 spans classified with the CGI, 29 spans, 
or 46%, were in the zone of the same classification 


Table 8. Number of spans classified with the CGI (proposed weights) by scenarios, and their adherence to 
the Laubscher’s adapted diagram. Note that scenario 4 is the most adherent. 

















Adherence to Laubscher’s Diagram (total 63 spans) Scenario 1 Scenario 2 Scenario 3 Scenario 4 
No. of adherent spans 17 24 26 29 
% of adherent spans 27% 38% 41% 46% 

















when inserted into Laubscher’s diagram (Table 8). 

With the best adherence verified for scenario 4, 
these proposed weights were considered theoreti- 
cal weights that will be used in the continuation of 
the development of the CGI. The 46% adherence 
between the CGI and the adapted Laubscher was 
considered extremely satisfactory because both 
have the same stability classification function and 
use the Bieniawski (1989) RMR index, but include 
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other different specific and peculiar variables and, 
therefore, higher adherence should not be expected. 
Finally, Table 9 summarizes how the theoretical 
weights were after the analysis, aiming the next 
step of field application in the mapping of the 27 
study caves, performing the last process in the CGI 
development flow to define the calibrated weights, 
those that will compose the official formulation of 
the index. 
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Table 9. Theoretical weight values for the classes of the four CGI variables. 






























































CGI 
CGI Variable Class Class Range 
(Theoretical Weights) 

Very poor 0 20 0 
; Poor 21 40 21 

Rock Mass Rating 
Regular 4] 60 39 

R 

Co Good 61 80 54 
Very Good 81 100 60 
e ’ Small Om 0.91 m 25 

Hydraulic Radius 
Regular 0.92 m 1.82 m 15 

HR 

om Large 1.83 m 3m 0 
a : Small 0m 3.31m 0 

Ceiling Thickness 
Cr Regular 3.32 m 7.64 m 3 
cy Large 7.65 m 10m ) 
ie Arch 10 

Ceiling Shape 
Planar 6 
(CS) 
Inverted Arch 0 
Legend: Best Class Regular Class Worst Class 





Calibrated Weights, Adjustment after Direct 
Field Application 

Adjustment after direct field application is the last 
process of the CGI development flow, when it goes 
for direct field application in cave mapping and 
span classification. The goal is to check whether the 
theoretical weights defined up to this stage are con- 
sistent with field reality or need to be adjusted. 


Mapping performed in the 63 spans of the 27 study 
caves showed that there was consistency in the 
CGI ratings for most caves, however, it was ob- 
served that in 26 spans of 17 caves, the CGI should 
have been stricter because structural fragility was 
observed, but did not reflect the index rating, that is, 
the ratings were better than they should be. 


Due to the need to increase the rigor of the CGI rat- 
ing, several reductions in the weights of all classes 
were tested until the new classification results were 
more compatible with the reality of field mapping. 


Table 10 presents the best alternative implemented 
for the calibrated weights, with a small reduction 
(between 15 to 30%) in the values of the RMR, 
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Ceiling Thickness and Ceiling Shape variable val- 
ues. The Hydraulic Radius variable did not change 
its weights because the tested reductions did not 
contribute to improvements in the final CGI rating 
and, therefore, this variable was not changed. 


After defining the calibrated weights, all 27 caves 
were again mapped in the field and classified by the 
CGI System demonstrating greater assertiveness 
between the classification results and the reality. 
Table 11 lists the 17 caves whose initial span CGI 
susceptibility ratings had not been adequate, con- 
firming the increase in the CGI rigor level, taking 
all ratings to a worse level than previously rated, 
including two really critical spans from the point of 
view of structural fragility mapped out from a high 
to very high susceptibility to structural instability 
rating. 


In summary, using the calibrated definitive weights, 
of the 27 caves in the study, 63% increased the in- 
dex sensitivity, and of the 63 spans, 41% increased 
their susceptibility to structural instability, better 


aligned with the reality observed in the field map- 
ping (Fig. 10). 
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Table 10. Reduction in the theoretical weight class values for the RMR, Ceiling Thickness and Ceiling 
Shape variables, defining the CGI calibrated weight (yellow highlighted lines). These new values are better 
adjusted with the mapping in the study caves. 






























































: CGI CGI Variation 
Cetvanable Cass (Theoretical Weight) (Calibrated Weight) (%) 
Very poor 0 0 
Poor 21 15 29 
Rock Mass Rating (RMR) Regular 39 30 23 
Good 54 45 17 
Very Good 60 60 
Small 25) DS) 
Hydraulic Radius 
Regular 15 15 
(HR) : 
Large 0 0 
Small 0 0 
Ceiling Thickness 
(CT) Regular 3 2 33 
Large 5) 5) 
Arch 10 10 
Ceiling Shape 
(Cs) Planar 6 4 33 
Inverted Arch 0 0 
Legend: Best Class Regular Class Worst Class 








27 Caves 63 Spans 





« CGI - Caves without susceptibility modification ® CGI - Spans without susceptibility modification 
® CGI - Caves with increased susceptibility to instability ® CGI - Spans with increased susceptibility to instability 


Figure 10. Comparative graphs after the remap of the 27 study caves using the calibrated weights. Note that from the 
27 caves, 63% had increased the rigor level, and in the 63 spans, 41% increased their susceptibility level to structural 


instability, more consistent with field reality. 
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Table 11. List of the 26 spans of the 17 caves classified with theoretical weights that were not consistent 


with reality, and the classification with the new calibrated weights that produced a more rigorous and closer 


to reality result. 













































































































































Results 

Example of Application of the CGI in a Cave 
Table 12 shows the CGI weights table for the class 
ranges of the 4 variables that are collected in the 


Cave. 


A digital spreadsheet was prepared to work on tab- 
lets or cell phones, facilitating data input and au- 
tomating the CGI calculations and classification 
during fieldwork. The study cave in this example is 
N4E_0022 which has 6 spans. Table 13 shows part 
of the spreadsheet, with the numbers in blue being 
the values of the classes of the 4 variables filled in 
during cave mapping. The respective CGI weights 


sy 


CGI classification CGI classification 
No. Cave Span 
(theoretical weight) (calibrated weight) 
1 N4E 0016 1 Low Moderate 
N4E 0022 if Low Moderate 
N4E 0022 2 Moderate High 
2 N4E 0022 3 Very low Low 
N4E 0022 4 High 
N4E 0022 6 Moderate High 
3 N4E 0023 1 Moderate High 
A N4E 0026 2 Low Moderate 
N4E 0026 6 High 
5 N4E 0027 1 Very low Low 
N4E 0037 1 Very low Low 
N4E 0039 3 Very low Low 
N4E 0044 1 Low Moderate 
: N4E 0044 2 Low Moderate 
9 N4WS 0008 3 Low Moderate 
N4WS 0008 5 Moderate High 
10 N4WS 0022 1 Very low Low 
11 N4WS_0050/51 5 Very low Low 
12 N4WS 0060 1 Low Moderate 
13 N5S 0001 3 Very low Low 
i N5S 0023 1 Very low Low 
NSS 0023 2 Very low Low 
il) NSS 0071 Il Low Moderate 
16 NSS 0073 1 Moderate High 
“a NSS 0075 1 Moderate High 
NSS 0075 D} Low Moderate 








are automatically inserted with the numbers in red. 
At the end, the sum of the weights leads to the CGI 
classification that indicates the level of susceptibil- 
ity to structural instability of the span as shown in 
Table 2. 


Finally, the map of the classified cave (N4E_0022) 
has its spans coloured with the standardized colours 
of Table 1 according to the CGI classification (Fig- 
ure 11). 


Accuracy Rate 
To measure the accuracy rate of the index, all 27 
study caves with their 63 spans were classified 
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Table 12. Class ranges of the 4 CGI variables with their respective weights 















































CGI Variable Class Class range Weight 
Very Poor 0-20 0 
Poor 21 - 40 15 
Rock Mass Rating (RMR) Regular 41 - 60 30 
Good 61 - 80 45 
Very Good 81 - 100 60 
Large 1.83 - 3.0m 0 
Hydraulic Radius (HR) Regular 0.92 - 1.82 m 15 
Small 0.00 - 0.91 m 25 
Inverted Arch 0 
Ceiling Shape (CS) Planar 4 
Arch 10 
Small 0.00 - 3.31m 0 
Ceiling Thickness (CT) Regular 3.32 - 7.64 m 
Large 7.65 - 10m 5 




















Table 13. Example of part of the CGI digital spreadsheet for the 6 spans of Cave N4E_0022. Note the values 
mapped in the field of the 4 variables (numbers in blue) and their respective weights and automatic CGI 
classification ( ). 












































CGI 
Suscepti- 
Ceiling Hidraulic 
Ceiling bility to 
Cave Thickness Radius 
RMR Shape Structural 
Span (CT) (HR) 
(CS) Instability 
(im) (m) 
of Cave 
Span 
51 
1 57 6.4 1.33 Planar 
Moderate 
36 
2 50 6.16 Zoe Planar 
High 
72 
3) 61 4.68 1.04 Arch 
Low 
Inverted 20 
4 37 7.97 2.35 
Arch Very High 
45 
5 34 8.4 1.11 Arch 
Moderate 
Inverted 32 
6 48 5.64 2.42 
Arch High 



































with the CGI and had continuous field monitoring 
throughout the research between 2015 and 2018 
during the mining operations. The field monitor- 
ing of the caves was performed by geostructur- 
al monitoring according to the cave breakdown 
mechanisms as described by Barbosa et al., 2019, 
from the first signs of damage until the total col- 
lapse when it occurred. It was planned to maintain 
the frequency of measurements always after mine 


Brandi: Cave Geomechanical Index (CGI) Geoconservation Research 


blasting. 

During the mine advancement operations towards 
the caves, the breakdown mechanisms were recog- 
nized and quantified as physical damage. The first 
results are shown in Fig. 12: 


= 8 caves (30%) showed no physical damage 
and remained stable and under monitoring; 
= 11 caves (40%) had direct collapse, that is, 


CGI=32 


Susceptibility to 
Structural 
CGI 
Instability of Cave 
Span 


CGI > 80 





Figure 11. CGI map of the study cave N4E_0022. 


showed no physical damage or any compro- 
mise of stability until the collapse, when op- 
erations were either very close to or above the 
cave; 

= 8 caves (30%) showed physical damage, 
2 of them remained stable and under moni- 
toring and six collapsed over time with the 
mining approach. 


iy 


In the analysis of the graph in Fig. 12, it is notewor- 
thy that 11 caves (40%) collapsed without the occa- 
sional physical damage common during mine ad- 
vancement, but only when blasting operations were 
on top of the cave, which confirms that iron caves 
are generally very resistant. If this amount is added 
to the 30% of caves that did not show any physical 
damage, 70% of the caves were found to be highly 
resistant and stable, although they had mining oper- 
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No Physical Damage (stable and under monitoring) 


With Direct Collapse 
With Physical Damage (stable and under monitoring) 


With Physical Damage (with later collapse) 


Figure 12. List of the 27 study caves and their physical damage status. Note that only 8 


caves registered physical damage (2 remaining stable and 6 collapsed over time) 


ations in their surroundings. 


To estimate the accuracy rate of the CGI, it was ob- 
served that the 11 caves with direct collapse had 16 
spans and did not show any physical damage until 
collapse. Of this total, only 1 span had been classi- 
fied as highly susceptibility to structural instability, 
but 15 were correctly rated as having moderate to 
very low susceptibility, that is, the accuracy of the 
index was 94%, indicating just during mapping, 
the low possibility of physical damage occurring in 
those respective spans (Fig. 13). 


On the other hand, considering the total of 63 study 
spans, the results were as follows (Fig. 14): 
= 27 spans (43%) showed no physical dam- 
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age and remained stable during monitoring; 
= 16 spans (25%) had direct collapse; 
= 20 spans (32%) showed physical damage. 


To estimate the accuracy rate, the only 8 caves with 
physical damage have a total of 26 spans, which 
only 20 were effectively damaged. The amount of 
the physical damage reached 68 records, that is, it 
was 68 times that this damage occurred repeated- 
ly inside the spans. Of this total, 76% occurred in 
spans correctly rated with CGI from high to very 
high susceptibility to structural instability, attesting 
to the reliability of the index (Fig. 15). 
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(calibrated weights) 


el eee N° of Caves with Direct Collapse 








[Nae occ? | 1 11 


| N4E 0039] 1 |  ~Low 
-Nae-00e | 2 N° of Spans of 11 Caves with Direct Collapse 


| N4E_0039 | 3 | Lowi 16 
| N4E_0061 | 4 | 











TN4E_oos6 | 1 | Me , | N° of Spans with Direct Collapse with CGI Susceptibility - Very Low/Low/Moderate 
15 194% 

N° of Spans with Direct Collapse CGI Susceptibility — High/Very High 
1/ 6% 











Figure 13. List of the 16 spans of the 11 caves that had direct collapse without showing physical damage, and the 
pyramid graph highlighting the accuracy rate of the CGI of 94%, where the spans were correctly classified with a 
Very Low / Low / Moderate susceptibility to structural instability. 
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Figure 14. List of the 63 spans and their physical damage status. 
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Caves with Physical Damage 
8 
N° of Spans of the 8 Caves 

26 

N° of Spans with Physical Damage 
20 

N° of Records of Physical Damages in Spans with Physical Damages 
68 
N° of Records of Physical Damages in Spans with CGI - Very Low/Low/Moderate 
16/ 24% 
N° of Records of Physical Damages in Spans with CGI - High/Very High 


52/ 76% 


Figure 15. Graph showing that of the 68 records of physical damage monitored in 20 spans, 76% occurred in spans 
correctly rated with a CGI of high to very high susceptibility to structural instability. 
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Figure 16. CGI maps of the only 8 caves that showed physical damage and their locations. 


























The accuracy rate of the CGI was 76% when com- 


The CGI was tested by classifying 63 spans of 27 
caves near mining operations in Carajas. All test 
caves were systematically monitored over 4 years, 
between 2015 and 2018, from the first signs of dam- 
age until their total collapse. 
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paring the locations of physical damage inside the 
spans, with their respective CGI classifications from 
high to very high susceptibility to structural insta- 
bility. On the other hand, the accuracy rate rises to 
94% for spans that were classified as moderate to 
very low susceptibility to structural instability and 


yi 





really did not show any physical damage. 


It is noteworthy that 11 caves (40%) collapsed with- 
out physical damage during the approach of the 
mining, and only failed when mining operations 
were on top of the cave, which confirms that iron 
caves are generally very resistant. If this amount is 
added to the 30% of caves that did not show any 
physical damage, 70% of the caves were found to be 
highly resistant and stable, although they had min- 
ing operations in their surroundings. 


Due to its pioneering application, focused only on 
iron caves, the CGI needs further trials and testing. 
It is possible that its variables and weights might be 
modified in the future. However, the first results are 
encouraging and motivate the authors to continue 
the evolution of the CGI and are already contem- 
plating the application of the CGI to karst environ- 
ments. 


Finally, the CGI outcomes open a path for safer 
mining operations close to the caves and can sup- 
port mine engineering, and the use of appropriate 
explosive charges, allowing the sustainability of 
mining and compliance with legislation in the spe- 
leological preservation. 
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